Abstract 



We study the Schwinger-Dyson equation for the fermion self-energy in 
massless and massive QED2, in the ladder approximation. When the fermion 
is massless (and the photon massless or massive), we check the reliability of 
this approximation by comparing its solutions with the exact ones. They agree 
only when the photon is massless. For a massive fermion and massless photon, 
we show that there is no consistent solution at all, the infrared divergences 
- introduced by the approximation forbidding even the trivial solution. When 

Q>^ . both fermion and photon are massive, we find a non-perturbative (extra) 

0\ I fermion mass generation (which survives in the limit when the bare mass of the 

fermion tends to zero). We argue that, in this case, the ladder approximation 
1-^ ■ will provide reliable results if the bare masses of the fields are large compared 

I't ^ , with the (dimensionful) coupling constant. 
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1 Introduction 



The Schwinger-Dyson Equations (SDE) of a model encode its perturbative and 
non-perturbative dynamics. Under some suitable approximation schemes, it is some- 
times possible to reduce them to a more manageable set of self-consistent equations 
involving the Green's functions of interest for the phenomenon under consideration. 
This approach has been extensively used for the study of non-perturbative effects 
in theories involving fermions and vector fields, with particular interest in the cal- 
culation of the fermionic self-energy, because of its relevance for (non-perturbative) 
chiral symmetry breaking and mass generation 0]. 

The simplest approximation for the calculation of the fermion self-energy is the 
so-called ladder approximation. It corresponds to summing over all the self-energy 
diagrams without fermion loops and without crossing between the different internal 
photon propagators. Although this is the simplest approximation, one can seldom 
extract any exact analytical result in this framework^. 

The different approximations used to simplify the SDE (large N, ladder) are not 
always easy to justify quantitatively, because of the lack of other calculations of the 
same objects by independent methods (except when some reliable lattice computa- 
tion is available). There are reasons to suspect that in low dimensions the popular 
truncations of the SDE could be in jeopardy. Indeed, the ladder approximation is 
essentially the Hartree-Fock method [|5[, which gives reliable results only when the 
fluctuations (i.e., current-current correlations) are negligible, and this could hardly 
be true in low dimensions. In particular, this approximation could predict a phase 
transition where there is none, as it happens when the mean-field approximation is 

used in the one-dimensional Ising model, for example. On the other hand, for some 

""^Considerable effort has been spent in applying this approximation to QED4. Maskawa and 
Nakajima [0 have proved that it predicts a self-consistent solution with a massive fermion. Fukuda 
and Kugo |j] converted the self-consistent integral equation to a differential equation plus boundary 
conditions, simplifying the analysis. 
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two-dimensional models one can solve the SDE exactly 0, so the possibility presents 
itself of comparing those solutions with the ones given by the ladder approximation. 

In this paper, we study analytically and numerically the SDE for the self-energy 
of the fermion in QED2 in the ladder approximation, exploring four different situ- 
ations regarding the bare masses of the fields: 

1) Both fields massless (Schwinger Model), 

2) Massless fermion and massive photon, 

3) Massive fermion and massless photon, 

4) Both fields massive . 

For the cases 1) and 2), we take advantage of the fact that the models are exactly 
solvable, which enables us to give a precise test of the approximation. We find that 
for 1) the approximation predicts the fermion mass to be zero, which is consistent 
with the exact solution 0. For 2) the numerical solution of the approximation 
predicts a non-zero mass, which is wrong (see Appendix). 

In the case 3), we show that there is no consistent solution at all, due to the wild 
infrared behaviour of the photon propagator. 

For 4), we present the numerical solutions of the corresponding equation, using 
the iteration method. We show that, when the bare masses of both fields are large 
compared with the coupling constant, the result provided by the ladder approxima- 
tion is reliable. 

The structure of the papers is as follows: in section 2 we introduce the model 
and obtain its corresponding ladder SDE for the fermion propagator. In section 3 we 
present our analysis of the four different cases, showing the corresponding numerical 

^Note that in the Schwinger Model there is chiral symmetry breaking, which manifests itself by 
a non-zero value of the chiral condensate {il^ip), but the fermion remains massless. 
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solutions, when they are non-trivial. In the Appendix we present the exact solution 
of cases 1) and 2). 



2 The model and its SDE 

All the cases we will consider can be obtained from the general Lagrangian 

£ = -\f^.F^^ + ^/i^A^A^ - \\{d^A^f + ^ - e 4 - m)V^, (1) 

where /i and m are the photon and fermion masses, respectively, e is the coupling 
constant, which has dimensions of mass. A is a covariant gauge fixing parameter 
when /i = 0. When /i 7^ 0, it improves the ultraviolet behaviour of the photon 
propagator without changing the physical quantities (the vector current is conserved 
quantum mechanically). It also allows us to take the massless photon limit safely. 

/^From eq.dH), we obtain the free fermion propagator 

S,{p) = , (2) 



and the free photon propagator 

k^ — /x^ k^ — fil 



G'-ik) = -^ C + ' (3) 



where we have defined fi\ = 

The SDE for the fermion self-energy is 

s-\p) = s,\p)-j:{p), (4) 

where S is the full fermion propagator, which we parametrize as 

^'f' = A(p) /- B(p) 
and E is the (IPI) self-energy. On the other hand, S is given by 

= {-^er J-0^,r S{k) r (p, k) G'^Xp - k) , (6) 
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where is the full vertex function and G'^^, is the full photon propagator. In the 
ladder approximation, is replaced by the tree vertex and G'^^ by the free 
propagator (^. Thus eq.(||) is approximated by 



■ 2 / d k 



m = i-teY J sik) Y g,Up - k) . (7) 

The set of eqs.(5) and (|^), when expanded in powers of e, corresponds to the series 
of 'rainbow' Feynman diagrams. They are equivalent to a set of two coupled integral 
equations for the scalar functions A and B: 

B = m-\ ^ — / dk—— , (8) 

2n A Jo A^k^ + 52 ^(p2 + ^2 + ^2)2 _ 4p2^2 

and 

_ e'^ 1 f°° kA ^(p2 + fc2)2 +^2(^p2 _^ _ 4p2^2 

-/I — 1 — ~ ~ I (lrC~ 



(9) 



An p'^fi'^ Jo A^k^ + 52 ^(p2 + ^2 + ^2)2 _ 4p2^2 

1 + A;2)2 + ^2^p2 ^ _ 4p2^2 



^ ^(p2 + A;2 + ^^2)2 _ 4p2^2 

(after Wick rotation and angular integration). Note that in the Feynman gauge 
(A = 1) eq. (j^) becomes just A = 1 and eq.(|^) yields 

B(p^)=m + — dk ^ ' , , , 10 

71 Jo P + ^(p2 + p + ^2)2 _ 4p2^2 

which is the only equation we have to consider. For latter convenience, we rewrite 
it in two different ways. First, defining x = and y = fc^, eq.(p!OD becomes 

B(x) = m + — [ dx^= ^ = ^^^Lx • (11) 

27cJ ./(x + ?/ + /x2)2_4a;^ ?/ + 52(/c2) 



On the other hand, for /i 7^ 0, we can simplify the treatment further by working 
in terms of a dimensionless function b such that B{p'^) = ^hijp'/ ^'^). In terms of 6, 
eq.([Tl|) becomes 

, / X r°° , 1 b(t) , , 

Us) = — + / dt , \i— . 12 



Now we study eqs. ([TI|) and ( |T2|) for the cases l)-4). 



3 Self-consistent solutions 



3.1 The case m = and /i = (Schwinger Model) 

In this case eq.(|lT]) reduces to 

B(.) = |l/<,, 1 (13) 

Although the kernel is singular at y = x, it might still be possible to get a finite 
answer for the rhs of eg. (pTS]) . Let us see that this is not he case. 

Assume that B{xo) is finite and non-zero in a neighbourhood of at least one 
point xq, then the Ihs of eq.(|T3|) at x = xq is, of course, finite and non-zero. But 
evaluating the rhs of (13) for x = Xq we get an infinite result, thus implying that 



B{xo) should be infinite, contradicting the original assumption. We conclude that 
the only solution of ([TT| ) is the trivial one: -B = 00. 

3.2 The case m = and /i 7^ 

As /i 7^ 0, we work in terms of the dimensionless function b of eq.(|T2D, which reduces 
to 

Hs) = r dt , ^ , . (14) 



We could not find any analytical solution of eg. (p^ , however, assuming that the 
solutions of eq.(p!4|) are decreasing functions, the following inequalities hold: 

which inserted into eq. (|T^ yield a constraint for 6(0) 

< 6(0) < 1^ . (16) 



^ An infinite B is not a solution, because it produces a on the rhs of ([ll]), and oo on the Ihs. 
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We have solved eq.(|T^ numerically by the usual iteration method, finding that 
there are non-trivial solutions for any (non-zero) value of the parameters. They 



consistent solution for different values of the parameters wich appear in the form of 
a dimensionless ratio g = In terms of g, (|1^) reads: < 6(0) < ng. 

In Figure 2, we plot the value of the self-energy at zero-momentum, as a function 
of the dimensionless coupling g. The numerical data fit very well with a power-law 
behaviour: 6(0) oc g'^^^. 

3.3 The case m ^ and fi = 

In this case eq.(^l]) yields 



and we see that the same argument which proves the non-existence of non-trivial 
solutions for the case 1), applies also here. However, the difference is that now, not 
even B = (or any constant) is a solution of (|l3)- Because of the infrared diver- 
gence in the photon propagator, the ladder series is ill-defined, and thus provides no 
consistent solution. In an improved approximation, one should take into account the 
fact that the photon propagator is modified by the fermion-loop correction. This 
modification will be particularly important in this case, because it will tame the 
infrared divergences by giving the photon a mass. We conclude from this example 
that neglecting the corrections to the photon propagator is not realistic in this case. 

3.4 The case m ^ and /i 7^ 

This case corresponds to the so-called Massive Schwinger Model, which has many 
interesting properties 0, and is not exactly solvable. From our point of view, 
we want to know which is the effect of an explicit bare fermion mass on the self- 



fall in the range predicted by ([T6|). In Figure 1 we show the behaviour of the self- 




(17) 
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consistent solution of eq.(|T^). 

Now we deal with eq.(p!2l) as it stands. From a discussion similar to the one of 

case 2), we get the constraint 



m < 6(0) < m + 



e2 



In Figure 3 we show the behaviour of the numerical solution of eq.(12) for g = 1. 
For fixed m, the profile of b is similar to the one of case 2), except for a vertical shift 
in m. In Figure 4 we show a plot of 6(0) as a function of the bare fermion mass. 
Note the it is very approximately linear with m. Let us consider this point with 
more detail. The explicit mass appears just as an additive constant in the r.h.s in 
eg . ([l^) . Then, in the first iteration, which is equivalent to a one-loop calculation, 
m just shifts the result corresponding to m = by a constant m//i. Repeating 
the iteration, m appears non-linearly so it would be naive to expect the solution 
for m 7^ to be just the one for m = shifted by a constant. But the numerical 
evidence shows that this is the case when m//i is larger than one. The reason is 
that when m is large compared to e, we can regard the ladder series as an expansion 
in powers of e/m, the leading term being the one-loop result. Thus we expect the 
non-linear part to be suppressed. To get an idea of when we can regard m as 'large', 
we estimate the relative ratio between the explicit and the dynamical masses after 
one iteration. 6(0) can then be evaluated analytically: 

2 

TP P 

m = - + -f{c) , (19) 

where 

/(c) = ^Inc , (20) 
— 1 

and c is the initial value of 6. The function /(c) has a maximum value / = | 
at c = 1. Thus, compared with the explicit term ^, the dynamical term gives a 
contribution less than irg for any initial value of 6. Then, when ^ ^ ng, we expect 
a behaviour approximately linear with m. 



Let us see now if in this limiting case the ladder approximation is reliable. We will 
approximate the exact SDE for very large m. Our approach is to replace G' and F 
in d^) by their large-m approximations. By large-m we mean that the dimensionless 
parameter e/m is very small, as well as the ratios p/m in the IPI functions. Let us 
consider first the correction to the photon propagator. Considering e/m << 1, it is 
easy to see that the leading correction is the usual fermion loop. For very large m, 
the contribution of this diagram is just 

n/.. = -^9,-^ (21) 

where we have also assumed that m is large compared with the external momenta. 
Then the full propagator corrected by this function will be just like the free one, 
but with /i^ replaced by fi"^ + e^/vr. 

Now, regarding the vertex function, we easily see that the first correction is of 
order 3 in e/m, and so we neglect it in the leading approximation, replacing F by 
the zeroth order function, i.e., F^ ~ 7^. Thus, putting together the approximations 
for G' and F, the result is that we get the ladder equation, but with a shifted photon 
mass. Now we can argue that the ladder approximation will be reliable if e/m << 1 
and e/fi « 1. The first condition is to justify the expansion, and the second to 
neglect the dynamical photon mass (|21|). 

Note that the assumption k/m << 1 in the computation of the fermion loop 
does not produce any harm at large k, because in that regime the behaviour of the 
vector field propagator is dominated by the k"^ in the denominator. 

Thus, we can also argue that under these conditions, the usual perturbation 
expansion should be reliable, i.e., the main contribution to the fermion self energy 
should be given by the one-loop diagram, or one-iteration with the bare mass as 
initial self-energy. The self-consistent calculation will improve this result, giving a 
contribution which will be small, but non-analytic in the coupling constant. In Fig. 
5 we show a plot of the one-iteration calculation, against the self-consistent solution 
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for a case where m = el\phx. Already the plots agree within a 5 percent error. 

We conclude this papers suggesting that the analysis we have made of this model 
could be extended to more realistic ones, where one might be interested in knowing 
if the non-perturbative effects are negligible or not in the computation of the full 
(bare + dynamical) fermion mass. Note that this aim is independent of the usual 
spontaneous chiral symmetry breaking study, because of the explicit breaking due 
to the bare fermion mass. 

Appendix: The exact and quenched solutions 

The exact solution of the Schwinger model has been known since the pionering work 
by Schwinger §]. We present here a calculation of the exact fermion propagator 
when the gauge field in the Schwinger Model is massive; this is a straightforward 
generalization of the massless case, but we show it here for the sake of completeness. 
It is also easy to discuss the effect of making the quenched approximation in this 
framework. In the path integral framework, the complete fermion propagator is 
defined by: 

S{x -y) = {iP{x)^{y)) = ^ JvAVijV^ Hx)i^{y) exp(« / d'xC) (22) 

where C is the Lagrangian of eq.(|l|), and is a normalization factor. The integration 
over the fermionic fields yields: 

S{x-y) = y PA det(^ + zm) (s I (^ + im)"^ I y) 

X exp[^ J dM-\F' + ^^^'A' - ^{d ■ A)')] . (23) 

Both the determinant and the inverse of the Dirac operator p = P + ie ^ can be 
calculated exactly: 

det(^ + im) = exp{-i J ^x^F'^'^F^,) , (24) 
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S{x -y) = exp{-ie[d "^d ■ A{x) - ed "^d ■ A{y) 

+ ^,,{d-\^^d^A,{x) - d-^^^^d.AM]} ■ (25) 

Inserting the results eqs.(]2^) and p5| ) into eg. (p3D , and performing the (Gaussian) 
integration over A, we obtain (now in Euchdean spacetime): 

S{x -y) = S'o(x - y) e^iW-^i e-^2(o)+A'.(x-j;) ^ ^26) 
with the definitions: 

( \ 2 f ^ ^ 1 ik-x 

/^From eg. (|26|) , one easily sees that the fermion remains massless since there is no 
extra pole at 7^ 0. 

It is very easy in this approach to understand what happens in the 'guenched' 
approximation. This is tantamount of discarding the determinant in eg. (^4]) , replac- 
ing it by 1. In diagrammatic language, it is eguivalent to summing all the diagrams 
which contain one electron line and any possible correction due to photon lines. 
Eguivalently, one can add to the rainbow diagrams the ones where the photons' 
propagators cross. Note that this approximation is gauge- invariant, what is trivial 
from the fact that we neglect a gauge-invariant guantity. This is achieved in the 
diagrammatic language by the crossed diagrams, which correct the vertex in such a 
way that the Ward identity is satisfied. 

In this approximation, we get also a massless fermion, as can be seen from eg . (p6D , 
disregarding the determinant's contribution. This leads to 

S{x-y) = So{x-y) (28) 

in the gauge A = 1. 
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Figure captions 

Figure 1 The self-energy b{s) for five different values of g = ( case 2 ). 

Figure 2 The zero-momentum self-energy 6(0) for some values of g, fitted by 
6(0) = ( case 2 ). 

Figure 3 The self-energy b{s) for = 1 and m — { case 4 ) . 

Figure 4 The zero-momentum self-energy 6(0) for g^ = 1 as a function of the 
ratio between the fermion and boson bare masses. The linear plot is the one-loop 
result ( case 4 ). 

Figure 5 The result of iterating the SDE (12) (for g — 1) once, starting from 

the initial value 6 = constant = — against the self-consistent solution of the same 
equation ( case 4 ). 
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